Cooperative long range protein-protein dynamics in Purple Membrane 
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Very recently, interprotein motions in a carboxymyoglobin protein crystal were reported from a 
molecular dynamics simulation [Phys. Rev. Lett. 100, 138102 (2008)]. We present experimental 
evidence for a cooperative long range protein-protein interaction in purple membrane (PM). The 
dynamics was quantified by measuring the spectrum of the acoustic phonons in the 2d Bacteri- 
orhodopsin (BR) protein lattice using inelastic neutron scattering. The data were compared to an 
analytical model and the effective spring constant for the interaction between protein trimers was 
determined to be k = 53.49 N/m. The experimental results are in very good agreement to the 
computer simulations, which reported an interaction energy of 1 meV. 

PACS numbers: 87.15.km, 87.16.dj, 87.16.D-, 83.85.Hf 



The high protein concentration in biological mem- 
branes might lead to long-range protein-protein inter- 
actions, on which there have been speculations, already 
some time ago Recently, interprotein motions in a 
carboxymyoglobin protein crystal were reported from a 
molecular dynamics simulation 0, Motions in pro- 
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FIG. 1: (Color online), (a) Sketch of the triple-axis scatter- 
ing geometry, qi | is the in-plane component of the scattering 
vector Q. (b) Reflectivity curve measured at T — 30° C. The 
dotted line is a fit of Lorentzian peak profiles including a q~ 4 
term, (c) BR trimers are arranged on a hexagonal lattice of 
lattice constant a « 62 A. (d) The interaction between the 
protein trimers is depicted as springs with effective spring 
constant k. 



teins occur on various length and time scales [4|, [5| , and 
the functional behavior of membrane proteins is likely 
to depend on the lipid bilayer composition and physi- 
cal properties, such as hydrophobic thickness and elas- 
tic moduli. How the variety of inter- and intra-protein 
motions, occuring over different time and length scales, 
interact to result in a functioning biological system re- 
mains an open field for those working at the interface 
of physics and biology. The dynamical coupling between 
proteins, i.e., cooperative protein dynamics, might be im- 
portant for the understanding of macromolecular func- 
tion in a cellular context because it might lead to an 
effective inter-protein communication. Here, we report 
collective inter-protein excitations in a biological mem- 
brane, the purple membrane (PM). 

Purple Membrane (PM) occurs naturally in the form 
of a two-dimensional crystal, consisting of 75% (wt/wt) 
of a single protein, Bacteriorhodopsin (BR), that func- 
tions as a light-activated proton pump, and 25% various 
lipid species (mostly phospho- and glyco-lipids) [6J. BR 
is a proton transporting membrane protein, formed of 
seven trans-membrane alpha-helices arranged around the 
photosensitive retinal molecule. The protein in the lipid 
matrix is organized in trimers that form a highly ordered 
2d hexagonal lattice with lattice parameter a f» 62 A, 
as depicted in Fig. Q] (c). The structure of PM is well 
established by electron microscopy, neutron and x -ray 
diffraction experiments reviewed in @, 0, H, @, EH U | . 
Here we present the unprecedented determination of col- 
lective protein-protein dynamics, i.e., acoustic phonons, 
of the 2d protein lattice in PM using coherent inelastic 
neutron scattering. The results allowed to determine the 
effective coupling constant. 

The experiments were performed on the INI 2 cold- 
triple-axis spectrometer at the Institut Laue Langevin 
(ILL, Grenoble, France). IN12 turned out to be highly 
suited for elastic and inelastic investigations in oriented 
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biological samples because of its flexibility, good energy 
resolution and extremely low background. It allows the 
measurement of diffraction and inelastic scattering in the 
same run without changing the set-up, which is crucial 
to assign dynamical modes to structural properties and 
molecular components. IN12 was equipped with its vac- 
uum box to avoid air scattering at small scattering an- 
gles and with vertically focusing monochromator and an- 
alyzer to increase the neutron flux at the sample position. 
There was no horizontal focusing, but the beam was col- 
limated to 40' — monochromator — 30' — sample — 30' — 
analyzer — 60' — detector. All scans were done in W- 
configuration with fixed kf — 1.25 A" 1 resulting in a q 
resolution of Aq — 0.005 A -1 and an energy resolution 
of AHlu = 25 /zeV. 

Deuterated PM was produced and hydrated by H2O in 
order to suppress the contribution of the membrane wa- 
ter to the phonon spectrum in coherent inelastic neutron 
scattering experiments and to then obtain information 
on the collective membrane dynamics. Because of the 
minuteness of the coherent inelastically scattered signals, 
the preparation of appropriate samples and experimen- 
tal set-ups is challenging in this type of experiments. We 
used a completely deuterated PM to enhance the col- 
lective protein-protein excitations over other contribu- 
tions to the inelastic scattering cross section. 200 mgs of 
deuterated PM suspended in H 2 was centrifuged down 
and the obtained pellet spread onto a 40 x 30 mm alu- 
minum sample holder. This was then partially dried to 
0.5 g water per gram of membrane over silica gel in a 
desiccator. The PM patches naturally align along the 
surface of the sample holder as they dry. The silica gel 
was then replaced by water, and the sample left to hy- 
drate to a lamellar spacing of d z = 65 A at 303 K (30°C). 

A sketch of the scattering geometry is shown in 
Fig. QJa). The experiments were carried out on PM 
stacks, i.e., PM samples with a regular repeat distance 
d z . The mosaicity of the sample (the distribution of nor- 
mal vectors with respect to the substrate) was checked 
by rocking scans to about 17 degs. The sample was 
aligned by centering the first reflectivity Bragg peak at 
q z » 0.1 A" 1 . Fi gure Qlb) shows a reflectivity curve 
measured at T = 30° C. From the three well developed 
Bragg peaks, the lamellar spacing d z was determined 
to be d z = 65.1 A, corresponding to an average inter- 
membrane water layer of 16 A, since the thickness of a 
dry PM fragment is 49 A. For the inelastic scans, Q was 
placed in the plane of the membranes (q\\)- 

The in-plane diffraction pattern of the 2d hexagonal 
protein lattice was measured and is shown in Fig. [5] 
Although the normal vectors of the membranes are 
well aligned with respect to the substrate, the (x, y)- 
orientation of the membrane layers is statistical and the 
signal a superposition of the different domains (powder 
average). All reflections can be indexed by a hexago- 
nal unit cell with a lattice parameter of 61.78 ± 0.73 A. 
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FIG. 2: (Color online). Diffraction pattern of the 2d protein 
lattice from qu = to 0.75 A" 1 . The reflections can be in- 
dexed by a hexagonal unit cell with a lattice parameter of 
61.78 ±0.73 A. The dashed line is a fit of the theoretical peak 
pattern to the data. The inset shows the third Brillouin zone 
in magnification. The arrows mark the positions of inelastic 
scans at constant qu values. 



Correlations and motions in membranes are often well 
separated in reciprocal space because of the largely differ- 
ent length and time scales involved. The prominent dis- 
tances in PM, such as lipid-lipid and BR-BR monomers 
and trimcrs for instance lead to spatially well separated 
signals. The same holds for the different time scales in- 
volved from the picosecond (molecular reorientations) to 
the nano- or microsecond (membrane undulations, large 
protein motions). The use of oriented samples further 
allows to separate correlations in the plane of the mem- 
branes, and perpendicular to the bilayers. Dynamics be- 
tween different protein trimers is expected to be dom- 
inant where the 2d BR diffraction pattern is observed, 
i.e., m a q\\ range of about 0.1 A" 1 to 0.6 A -1 . Be- 
cause elastic and inelastic scattering at small momentum 
transfers was dominated by the rather strong (1,0) and 
(1,1) reflections, systematic inelastic scans were taken at 
q\ 1 -values in the third Brillouin zone of the 2d pattern, be- 
tween 0.34 A -1 and 0.46 A" 1 . The inset in Fig.[2]shows 
the Bragg peaks in the third Brillouin zone in magnifi- 
cation and marks the positions of inelastic scans at con- 
stant g|| values. Figure [3] shows a constant-^ scan at 
<7|| = 0.39 A -1 . The total inelastic signal typically con- 
sists of a Gaussian central peak due to instrumental res- 
olution, quasielastic broadening, which is described by a 
Lorentzian peak shape, and, for this qu values, two pairs 
of excitations. Because of the pronounced and symmet- 
ric inelastic signals, we are sure that the observed peaks 
are not spurious effects. The quasielastic signal most 
likely contains contributions from coherent and incoher- 
ent scattering, i.e., from auto- and pair-correlated molec- 
ular motions. 
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FIG. 3: (Color online). Energy scan at q\\ = 0.39 A -1 . The 
total inelastic signal consists of a central peak (Gaussian), 
a Lorentzian quasielastic contribution, and the excitations, 
described by damped harmonic oscillators. 



Propagating modes are experimentally difficult to ob- 
serve in the scattering function S(q, u>) as they often show 
up as a small peak or a soft shoulder in the tails of the 
central mode and quasielastic broadening, rather than as 
a clear excitation maximum. In contrast, in the longi- 
tudinal spectrum defined by Ci(q,u>) = (uj 2 / q 2 )S(q,Lo), 
the multiplication by u> 2 suppresses the central mode and 
quasielastic scattering and makes the excitations (prop- 
agating sound modes) more easily visible for analysis. 
Systematic inelastic scans for qu values between 0.34 A -1 
and 0.46 A" 1 are shown in Figure HJa) ■ The position of 
the Bragg reflections are also marked in the figure. 

In order to interpret the data, the excitation spectrum 
of the 2d protein lattice was modeled analytically. In 
a simple approach, the protein trimers were taken as 
the centers of a primitive hexagonal lattice with lattice 
constant a = 62 A, and the spectrum of the acous- 
tic phonons was calculated. The model is depicted in 
Fig. [Ijc). The basic hexagonal translations are marked 
by arrows. The interaction between the protein trimers 
is contained in springs with an effective (longitudinal) 
spring constant k (Fig. QJd)). The calculated Ci(q,u) is 
shown in Fig. HJb). The statistical average leads to a su- 
perposition of the different phonon branches, which start 
and end in the hexagonal Bragg peaks (at Hw = 0). The 
absolute phonon energies can not be determined from the 
model, but depend on the coupling constant k. So the 
energy of the phonon curves in Fig. [Jb) was scaled to 
match the experiment. Note that because the proteins 
trimers were treated as dots with an effective mass of 
M tr , the calculation does not include any contributions 
from intra-protein or intra-trimer dynamics, i.e., possible 
optical modes and phonons. 

To compare to the data, the most pronounced phonon 
branches from Figure |U[b) were plotted as solid lines in 




FIG. 4: (Color online), (a) Experimental Ci(q,uj), measured 
in the the third Brillouin zone, (b) Calculated excitation spec- 
trum Ci (q, ui) in the range of the experimental data. The most 
strongest phonon branches are also shown in part (a) as solid 
lines. 

part (a). In the qii range between 0.34 A -1 and about 
0.43 A -1 , the experiment well reproduces the calculated 
phonon curves. The agreement is less good for g| values 
above 0.43 A -1 . Only the strongest phonon branches 
are visible in the data, the weaker branches can most 
likely not be resolved from the background. The intra- 
protein dynamics can most likely be observed at higher 
q\\ values (because of the smaller distances involved), 
around the corresponding correlation peaks, and should 
in a first approximation not interfere with the inter pro- 
tein dynamics. So it is reasonable to expect that the 
main structure of the calculated branches is found in 
the data, but the calculation will not reproduce a pos- 
sible fine structure due to internal motions. The arrow 
in Fig. Ufa) marks the position of the energy scan in 
Fig. [3l at gn = 0.39 A" 1 . The scan cuts the theoret- 
ical dispersion curves twice, at energy values of about 
hui = 0.7 meV and hui = 0.9 meV. The data in Fig. [3] 
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were fitted using two damped harmonic oscillators at 
hu>o = ±0.71 meV and Tiloq = ±0.92 meV, in very good 
agreement to the calculations. It was difficult to observe 
excitations at energies below about 0.5 meV because of 
the limited instrumental resolution and the quasielastic 
scattering contribution. The excitations would then show 
up as shoulders of the central peak and quasielastic scat- 
tering, only, rather than as pronounced inelastic peaks. 

Experiment and calculation can be compared to the 
MD simulations in Ref. [2|, where a peak in the energy 
spectrum at 1 meV was identified as a translational inter- 
molecular proteimprotein interaction vibration in a car- 
boxymyoglobin protein crystal. The energy agrees well 
with the energy of the zone boundary phonon in PM of 
1.02 meV, as shown in Fig. [H but no g-dependence was 
determined in the computational work. The good agree- 
ment of the energy values in the two systems most likely 
stems from the very high protein density in PM, which 
makes it almost crystal like. The computational work 
thus strongly supports the interpretation of our data as 
collective proteimprotein excitations. 

The commonly assumed interaction mechanism be- 
tween inclusions in membranes is a lipid-mediated in- 
teraction due to local distortions of the lipid bilayer 
0) 13, 14, [ID, [l6j]> with a strong dependence on the 
bilayer properties, in particular elastic properties. The 
PM might, however, be a special case because there are 
very few lipids between neighboring BR proteins [l7j |. 
While the nature of the interaction still will be mainly 
elastic, it is not likely to be purely lipid-mediated but for 
the most part a direct protein-protein interaction. The 
strength of the interaction can be determined from the 
data in Fig. [4j The energy of the zone-boundary phonon 
at the M-point of the hexagonal Brillouin zone (for in- 
stance at a <7| | value of 0.35 A" 1 ) relates to the coupling 
constant by M tr uj 2 — 6k. Because this energy is deter- 
mined as Hlu — 1.02 meV, the effective protein-protein 
spring constant k is calculated to k = 53.49 N/m [2l|. 
On the microscopic level, displacing the BR trimer by 1 
A yields a force between neighboring trimers of 5.3 nN. 
There is therefore strong protein-protein communication 
in PM. Using the same approach, the spring constant for 
graphite for comparison is calculated to 27,000 N/m for 
the in-plane interaction, and 3.5 N/m for out-of-plane in- 
teractions. The force constant that we measure in PM 
thus is 1-2 orders of magnitude larger than the effec- 
tive van-der-Waals force constant in graphite, but 2-3 
orders of magnitude weaker than a C-C bond. The mech- 
anism might be relevant to model the photo cycle in PM 
where the BR proteins undergo small structural changes 
which then propagate to neighboring proteins by 
the coupling k. Because the energy of the excitation is 
about 1 meV w 12 K, this process is highly populated at 
room temperature (96%) and might therefore be highly 
relevant for biological function. 

In conclusion we present experimental evidence for 



a cooperative long range protein-protein interaction in 
purple membrane. The effective spring constant for the 
interaction between protein trimers can be determined 
from the acoustic phonon branches as k = 53.49 N/m. 
The finding might be relevant for protein function, i.e., 
for the photo cycle of BR proteins, where the BR 
monomers undergo small structural changes. Future ex- 
periments will address cooperative intra-protein-trimcr 
and also intra-monomcr dynamics. To make a clear re- 
lationship to protein function, protein dynamics of acti- 
vated proteins, i.e., proteins undergoing the photo cycle 
will be studied. 
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